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ABSTRACT 

Recent observations and theoretical investigations of neutron stars indicate that their atmo- 
spheres consist not of hydrogen or iron but possibly other elements such as helium. We cal- 
culate the ionization and dissociation equilibrium of helium in the conditions found in the 
atmospheres of magnetized neutron stars. For the first time this investigation includes the 
internal degrees of freedom of the helium molecule. We found that at the temperatures and 
densities of neutron star atmospheres the rotovibrational excitations of helium molecules are 
populated. Including these excitations increases the expected abundance of molecules by up to 
two orders of magnitude relative to calculations that ignore the internal states of the molecule; 
therefore, if the atmospheres of neutron stars indeed consist of helium, helium molecules and 
possibly polymers will make the bulk of the atmosphere and leave signatures on the observed 
spectra from neutron stars. We applied our calculation to nearby radio-quiet neutron stars with 
Bdipoie ~ 10 13 -10 14 G. If helium comprises their atmospheres, our study indicates that isolated 
neutron stars with T BB ~ 10 6 K such as RXJ0720.4-3125 and RXJ 1605. 3 +3249 will have He + 
ions predominantly, while isolated neutron stars with lower temperature (Tbb ~ 5 x 10 5 K) 
such as RXJ 1856.5-3754 and RXJ0420. 0-5022 will have some fraction of helium molecules. 
We found that ionization, dissociation and electric excitation energies of helium molecules 
are larger than 100 eV at B > 10 13 G. On the other hand, rotovibrational excitation energies 
are in the range of 10-100 eV at B — 10 12 -10 14 G. If helium molecules are abundant, their 
spectroscopic signatures may be detected in the optical, UV and X-ray band. 
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1 INTRODUCTION 



Hydrogen has been considered as the surface composition of 
isolated neutron stars (INSs) because gravitational stratifica- 
tion forces the lightest el ement to the top of the atmosphere 
dAlcock & Illarionovl [l980h . Only a tiny amount of material is 
required to co nstitute an optically thick l ayer on the s urface 
jRomanilll987h. However, r ecent studies of IChang et al. I l l2004l) 
and lChang & Bildstenl Jiooj) have shown that the NS surface may 
be composed of helium or heavier elements since hydrogen may 
be quickly depleted by diffusive nuclear burning. Observation- 
ally, helium and heavier element atmospheres have been proposed 
for interpreting the spectral features observed in several INS par- 
tially because the existing hyd rogen atmosphere models do not re- 
produce the observed spectra I Sanwal et alj|2002k lHailev & Moril 
l2002l : lvan Kerkwiik & Kaplanll2006T) . However, atomic and molec- 
ular data in the strong magnetic field regime are scarce for non- 
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hydrogenic elements. Accurate at omic and molecular data ar e 
available mostly for the He + ion l Pavlov & Bezchastnov 120051) . 
the helium atom dNeuhauser et alJl987l : lDemeur et al] 19941). He 3+ 
i Turbiner & Lopez Vievral 1 20041) and He 2+ I Turbiner & Guevara 
l200rj) . Helium molecular binding en ergies have been cru dely cal- 
culated by density functional theory jMedin &~La i 2006a. b) (here - 
after ML06). Unlike hydrogen atmospheres dLai & Salpeteill 19971) . 
the ionization and dissociation balance in strongly-magnetized he- 
lium atmosphere has not been investigated yet. 

In this paper, we ex tend our Hartree-Fock type calcula- 
tion l lMori & Haiievl 120021) to helium molecules in the Born- 
Oppenheimer approximation. For molecular ions that exist in 
strong magnetic fields (B = 10 12 — 10 14 G), we achieved S \% 
and .5 10% agreement in binding energies and vibrational energies 
in comparison with other more accurate studies mainly on hydro- 
gen molecules. Including numerous electronic, vibrational and ro- 
tational states, we studied ionization and dissociation equilibrium 
in helium atmospheres at B = 10 12 -10 14 G. We also applied our 
calculations to several INSs which may have helium atmospheres 
on their surfaces. 
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2 MOLECULAR BINDING AND VIBRATIONAL 
ENERGY 

At first we adopt the Born-Oppenheimer approximation and ne- 
glected any effects associated with motion of atoms and molecules 
in a magnetic field. Later, we will discuss rotovibrational states (^3) 
and how the finite nuclear mass modifies results (§5$. In the Lan- 
dau regime (J3 Z > 1 where B z = B/(4J x 10 9 Z 2 G) and Z is the 
atomic number), bound electrons in an atom and molecule are well 
specified by two quantum numbers (m, v). m is the absolute value 
of a magnetic quantum number (which is negative to lower the total 
energy in strong magnetic fields) and v is a longitudinal quantum 
number along the field line. We consider only tightly-bound states 
with v = 0. Electronic excited states with v > have small binding 
energies, therefore their population in the atmosphere is tiny due 
to small Boltzmann factors and pressure ionization. Hereafter we 
denote atomic and molecular energy states as (/?;i,m 2 ,m3,m4, ...). 

We computed molecular binding energies with a simple 
modification to our Hartree-Fock type calculation for atoms 
dMori & Haiievl l2002h . We replaced the nuclear C oulomb term 
V m (z) in the Schrodinger equation by dLaietalJl992h 



V m (z) ~- 
where 

VJz) 



(1) 



(2) 



The function O,,, is the ground state Landau wavefunction, and a is 
the separation between two nuclei. We added as the Coulomb 
repulsion energy between the two nuclei. We computed binding 
energies with a grid size Aa ~ 0. 1 [a.u.] up to a ~ 1 [a.u.] and 
Aa ~ 0.01 [a.u.] near the energy minimum. Figure [T] shows the 
binding energy cur ve of He2 at B = 10 12 G fitted with the Morse 
function defined as jMorsell 19291) 



U(a) = D m {\ - exp [-((a - a a )]} 2 - E„, 



(3) 



where £,„(> 0) is the molecular binding energy for an electronic 
state m (in this paper this usually denotes the magnetic quantum 
number of the outermost electron), a is the separation between two 
nuclei at the minimum energy. We defined two different dissocia- 
tion energies: D,„ = E,„ - 2E a (gs) and £>„, = E„, - E a (m a ,mi,) - 
E a {m c ,m d ). E a (gs) is the ground state energy of an atom (e.g., 
(0, 1) state for Helium atom) and E a (m a , nib) is the energy of an 
atom in the (m„, mt) state. Each of the atomic m quantum numbers 
{m a ,m b ,m c and m d ) corresponds to one of the molecular m quan- 
tum numbers (wi,W2,W3 and 0J4) so that E a {m a ,m\,) + E a (m c ,ma) 
is the smallest. For instance, helium atoms in (0, 3) and (1,2) state 
are the least bound system into which He 2 in the ground state (i.e. 
(0, 1, 2, 3) state) will dissociate. Note that a molecule dissociates to 
atoms and ions when E„, < 2E a (gs), while the molecular binding 
energy approaches E a (in a ,m b ) + E a {m c , m d ) at large a. 

The calculated binding energy values are not smooth near the 
energy minimum (Figure [TJ. This is due to our numerical errors. 
Binding energy does not change by more than 0.1% for Aa = 0.01 
[a.u.] near the energy minimum. We determined E„, from the fit- 
ting procedure using the function given by equation $3$ since we 
found that it provides more accurate results than the minimum en- 
ergies from our grid calculation. However, in most cases, E,„ from 
our grid calculation and the fitted E„, do not differ by more than 
1%. We computed D„, and D„, using the atomic data we calculated 
numerically. 

Our results for Hi and H2 are in good agreement 
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Figure 1. Binding energy curve of Heo at B = 10 12 G. The solid curve 
is the fitted Morse function. The dashed line corresponds to the energy of 
two Helium atoms in the ground state (0, 1) (= -E,„ + D,„). The dotted line 
corresponds to the summed energy of one Helium atom in (0, 3) state and 
the other in (1, 2) state (= ~E m + D m ). 



Table 1. Total binding energy [eV] of Hj (left) and H 2 (right). 
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Notes : LS96: lLai & Salpetei] jl996h . TL03: iTurbiner & Lopez VievrJ 
120031) 



with ITurbiner & Lopez Vievr j d2003h (hereafter TL03) and 
lLai & Salpetej jl996h (LS96) with less than 1% deviation in to- 
tal binding energy (Tab. [T}. TL03 performed highly accurate vari- 
ational studies mainly on one-electron molecular systems (e.g., 
LS96 studied hydrogen molecular structure simi- 



Ht, H 



He3+). 



larly by a Hartree-Fock calculation in the adiabatic approximation. 
While our calculation takes into account higher Landau levels using 
perturbation theory, the difference in binding energies by includ- 
ing higher Landau levels is tiny for helium atoms and molecules at 
B ^ 10 12 G. 

Similar to hydrogen, the ground state configuration is 
{m\,m2) = (0,1) for He 2 / and (mi, m2,m^,mn) = (0,1,2,3) for 
He 2 at B ^ 10 12 G. The accurate comparisons for helium molecules 
are Demeur et al. Ul994h and ML06 who computed He 2 binding en- 
ergy by Hartree-Fock theory (table|2]l. Our results agree with ML06 
within 1%. ML06 also computed helium molecular binding ener- 
gies using density functional theory (ML06)Q However, their DFT 
results are less accurate than those of Hartree-Fock calculation; the 
binding energies are overestimated by ~ 10% (ML06). 



1 Although the results of ML06 are mostly from density functional calcu- 
lation, they showed s ome results from Hartree-Fock calculation based on 



Lai & Salpetej 0996) for comparison. 
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Table 2. Total binding energy [eV] of He2 in comparison with ML06. 
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Figure 2. Binding energy curve of He2 in (0, 1,2, m) states at B = 10 12 
G. The asterisk points are the binding energies from our numerical calcu- 
lations. The dashed line corresponds to the energy of two Helium atoms in 
the ground state. 



perturbation theory, [khersonskii i ll 985b has shown that the cou- 
pling energy is tiny (less than 1%) compared to the total binding 
energy. Neglecting the coupling, the rotovibrational energy levels 
are well approximated by 



£nav - £v + Cjva- 

ev is the aligned vibrational energy given by 



2 4D m 



1 

V+ 2 



(4) 



(5) 



The integer 0) is the quantum number for the aligned vibra- 
tion and Ho)» is the aligned vibrational energy quanta jMorse|[T929l : 
iKhersonskil 19851) . 

On the other hand, the transverse rotovibration energy (6jva) 
consists of transverse vibr ation and rotation around the magnetic 
field axis and it is given as dKhersonskiilll985l) 



6/VA = fko, \N + 



IM + 1 



(6) 



The integer N(^z 0) is the quantum number for the transverse vi- 
bration, while the integ er A is the projectio n of angular momentum 
in the B-field direction dKhersonskiilll985l) . hQ. B is the nuclear cy- 
clotron energy (= ZheB/(Am p c) = 63(Z/A)B l2 [eV] where A is 

the atomic mass) and a>, = ^4ll> 2 ± + Q? B . a> ± is the transverse vibra- 
tional energy quanta. The nuclear cyclotron energy term takes into 
accou nt the magnetic restoring force on the nuclei dLai & SalpeteJ 
Il996h . In the following subsections, we calculate vibrational energy 
quanta in the Born-Oppenheimer approximation. 



2.1 Electronic excitation 

The electronic excited states in molecules occupy higher (m, v) 
states than those in atoms. Since excitation energies from the 
ground state to (0, 1,2, m) states with m > 3 are small, there may 
be numerous tightly-bound electronic excited states until they dis- 
sociate into atoms and ions at large m. We did not consider excited 
states with v > because their binding energies are small therefore 
they are likely to be dissociated. 

We calculated binding energies for the (0, 1,2, m) state up to 
m = 9 and estimated binding energies for higher m states using the 
well-known m dependence of the energy spacing AE„, ~ m (^f) 
dLai et al ][l992). We found that the difference between the exact 
solutions and those from the scaling law is tiny at m > 9. Figure 
[2] shows the He 2 binding energy of (0, 1,2, m) states at B = 10 12 
G. Note that the excited states of He2, (0, 1, 2, m) with m > 5, are 
unbound with respect to two atoms in the ground state at B = 10 12 
G. 



3 ROTOVIBRATIONAL EXCITATION 

We consider molecular excitation levels associated with vibrational 
and rotational motion of molecules in a magnetic field. In con- 
trast to the field-free case, the strong magnetic field induces molec- 
ular oscillations with respect to the field line similar to a two- 
dimensional harmonic oscillator. Accordingly, there are three types 
of molecular motion; vibration along and transverse to the magnetic 
field and rotation around the magnetic field. Hereafter we briefly 
describe energy levels of rotovibrational states. 

Stri ctly speaking, the aligned and transver se vibrations are 
coupled jKhersonskiilll984lLai & Salpetenll996l) . However, using 



3.1 Aligned vibrational excitation 

In the Born-Oppenheimer approximation, the motion of two nuclei 
along the magnetic field is governed by the binding energy curves 
determined in section [2] It is therefore straightforward to calculate 
aligned vibrational energy quanta using the results from section [2] 
We fit the Morse function given by equation (3j to molecular bind- 
ing energies as a function of the nuclear separation a. Once f is de- 
termined, the aligned vibrational energy quanta is given by (LS96) 

fuo, = n\^A 2 (7) 

where ft is the reduced mass of the two nuclei in units of the 
electron mass (918 for H 2 and 3675 for He 2 ). At large a, another 
electr on configuration is mixed with tightly-bound states dLai et al .1 
1 19921) . Since configuration interaction is neglected in our calcula- 
tion, our £ values (therefore hu>\\) are overestimated by 10-30% in 
comparison with LS96 (table[3}. We found that { is nearly identical 
for different (tightly-bound) electronic excited states. Therefore we 
computed fra)\\ for electronic excited states with large m (for which 
we did not perform grid calculation) using f from the lower excited 
states. In most cases, our aligned vibrational energy quanta agree 
with other more accurate results w ithin ~ 10% (table [3 ) . Tabl e [4] 
compares our results for He^ with lTurbiner & Guevaral J2006h at 
various magnetic fields. Table [5] shows aligned vibrational energy 
quan ta for helium molecular ions alon g with some results on He, + 
from lTurbiner & Lopez Vievral d2004l) . 

It is apparent that the discrepancies with other calculations are 
larger for helium molecules than for hydrogen molecules. There are 
two effects both of which reduce the accuracy of our results particu- 
larly for highly-ionized molecules at low magnetic fields. First, we 
did not take into account configuration interaction, while the work 
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Table 3. Aligned vibrational energy quanta tuj^ [eV] of Hj (left) and H2 
(right). 
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Notes: We multiplied the results of TL03 by a factor of two to correct 
the discrepancy due to different definitions of aligned vibrational energy 
quanta. 



Tab le 4. Comparison of aligned vibrational energy quanta fou)\\ [eV] of He^ 
with lTurbiner & G uevara (20q3). 
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Notes: Similar to table[3] we multiplied the results of TG06 by a factor of 
two. 

by Turbiner et al. employed the full 2-dimensional variation en- 
ergy calculation. The degree of configuration interaction is larger 
at small f} z since the increasing effects of the nuclear Coulomb 
field mix different electron configurations. Second, highly-ionized 
molecular ions are either unbound or weakly-bound at low mag- 
netic fields, therefore the numerical errors in binding energies sig- 
nificantly affect determination of £ since the binding energy curve 
is shallow. Nevertheless, the accuracy of vibrational energies of 
highly-ionized molecules such as He^ is irrelevant for dissocia- 
tion balance since they are unbound or their abundance is negli- 
gible (section[6]l. For the abundant molecular ions in 10' 2 — 10 14 G 
(e.g., He} and He2), we expect the accuracy of aligned vibrational 
energy is < 10%. 

3.2 Transverse vibrational excitation 

We calculated the energy curve as a function of transverse position 
of nuclei R following Ansatz A described in section IIIB of LS96. 
We fixed a to the equilibrium separation a and supposed that the 
two nuclei are located at (±R/2, ±ao/2). As LS96 pointed out, this 
method is appropriate for small R (< p where p is the cyclotron 
radius) and gives only an upper limit to transverse vibrational en- 
Table 5. Aligned vibrational energy quanta fmja [eV] of helium molecular 
ions. 
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Notes: The numbers in t he brackets show the results from 
iTurbiner & Lopez Vievri] ]2004l) . 



Table 6. Transverse vibrational energy quanta hw± [eV] of (left) and 
H2 (right). The numbers in the brackets are transverse vibrational energy 
[eV] computed using the perturbation theory. 
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Table 7. Transverse vibrational energy quanta [eV] tuo± of helium molecu- 
lar ions. The numbers in the brackets are transverse vibrational energy [eV] 
computed using the perturbation theory. 
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(a) Transverse vibrational e nergy quanta of He^ ions from 
ITurbiner & Lopez Vievri] 120041) are 11 eV and 51 eV at B l2 = 1 
and 10 respectively. 

ergy quanta hco ± . We replaced the nuclear Coulomb term in the 
Schrodinger equation Viz) by 

^M>-K.(«-a.fM»?.!) 

where 

V m (z,R/2)= f c?A^^. (9) 
•J \?x - R/2\ 

We added , I .,, as the Coulomb repulsion energy between two 

(al+R 2 )" 2 * bJ 

nuclei. Once we calculated the molecular binding energy at dif- 
ferent R grid points, we fit a parabolic form ^poj^R 1 to binding 
energies at R < p. We found that fko± is nearly identical for differ- 
ent electronic excitation levels. Therefore, we adopted fia> ± of the 
ground state for electronic excited states. Our transverse vibrational 
energy quanta agree with those of LS96 and TL03 within 10% (ta- 
ble[6j. Table|7]shows the results for helium molecular ions i n com - 
parison with those for He^ + from lTurbiner & Lopez Vievral ( 120041) . 
Compared to aligned vibrational energy, our results for transverse 
vibrational energy are in better agreement with other more accu- 
rate results. The better agreement is well-understood because the 
transverse potential well is deeper than in the aligned direction in 
magnetic field, therefore our results are less subject to the numeri- 
cal errors as discussed in section [3T| 

3.2.1 Perturbative approach 

It is also possible to estimate the transverse vibrational energy per- 
turbatively, by calculating the lowest order perturbation to the en- 
ergy of the molecule induced by tilting it. We assume that the en- 
ergy of the tilted molecule is almost the same as the energy required 
to displacing the electron cloud relative to the molecule by an elec- 
tric field (E) : 

e[ l) = (K\eEx\K) = 0. (10) 
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The first-order change to the wavefunction is 



„(i) 



e 



(k' \eEx\ k) 



(11) 



where e£ 0) is the unperturbed energy of the state k and k' denotes 
the other states of the system. 

The expectation value of x for this situation is 

\{k'\x\k)\ 2 



r sin 9 = (x) 



- £ E 



,(0) _ (0) 



(12) 



where r is half the distance between the nuclei. We can solve for the 
value of eE that we need to apply to give a particular displacement 
and substitute it into the expression for the energy of the state to 
second order 

-i 



JZ) 



( ^) 2 Z 



K*'|xk)r 

JO) 



r sin" 8 



E 



\(k'\x\k)\ 



2 \ 



JO) 



(0) 



(13) 



The only states that contribute to the sum have m' = m ± 1 for the 
electronic states of the molecule (the rotational states do not count 
because that is what we are examining). 

States that have been excited along the field (v > 0) or by 
increasing the Landau number do not have much overlap in the 
integral in the numerator and also a large energy difference in the 
denominator. 

The frequency for low amplitude oscillations is given by 

J_(y \{k'\x\k)\ 2 ^ 



I 



E 



Kir* \x\k)\ 1 



JO) _ (0) 



J0) _ (0) 



(14) 



where / is the moment of inertia of the molecule, 2Mr 2 (M is the 
nuclear mass). The size of the molecule has cancelled out. 

For a single electron system in the ground state if we assume 
that the bulk of the contribution to the sum in the equation is given 
by the m = 1, v = state we have 

V2, 
P 



y2 /e m =i,v=o - e m =o,v=o \ 1/2 i, c ., >i 

J ± ~ —T- y ^ J |\/m=l,v=ol/m=0,v=0/| 



(15) 



where the final term is the overlap of the longitudinal wavefunction 
of the two states. 

For a multielectron system, evaluating equation d 1 3 1 > and d!4t 
is somewhat more complicated. For clarity of nomenclature we 
shall write the wavefunctions of the various electronic states of 
the molecule as %' . The symbol % denotes the state that we are 
focused upon. The change in the energy of the system due to an 
applied electric field is 



{ r K\NeEx\K) = Q 



where N is the number of electrons and 



(16) 



(17) 



where j counts over the electrons in the molecule. We have as- 
sumed that the multielectron wavefunctions are normalized such 
that {K\K) = 1 

The first-order change to the wavefunction is 

(<K' \NeEx\7Q 



= E 



M. 



,0, 
fc 'A- 



-<K' 



(18) 



'A" 

where e^l is the unperturbed energy of the state %' . 

Now let us calculate the expectation value of x for this situa- 
tion, 



(x) = NeE V 



■K' 



\(K' \x\K)\- 

(0) _ (0) 



(19) 



where r is half the distance between the nuclei. We can solve for 
the value of NeE that we need to apply to give a particular displace- 
ment and substitute it into the expression for the energy of the state 
to second order 



(NeE) 2 ^ 



\(K' \x\<K)\ 2 



E 



lOTC'i-wr 



(20) 

C 'A" C 'K 

In strongly magnetized atoms or molecules, it is natural to ex- 
pand the wavefunctions in terms of the ground Landau level us- 
ing various values of m. We assume that the wavefunction for each 
electron is written as 

In this case we have 



(21) 



(kMkj) = P - V2(m,- + 1) (/■'/}) if m, = mj + 1 



(22) 



otherwise it vanishes. 

Combining these results yields an estimate for the frequency 
of low amplitude oscillations of 



yvV2 



1/2 



(23) 



/ A , e /»+i - f »; \ 1/2 



(24) 



p \M(m + 1) j 

where the subscript on energy labels the value of m in the outermost 
shell. The number of electrons appears in this equation because the 
expectation value of x. is a factor N smaller than the expectation 
value of x for the single shifted electron. 

For the ground state of the molecule we have m + 1 = N yield- 
ing a simpler expression, 

V2, 
P 

Table [6] and [7] compare the perturbative estimates with the numeri- 
cal calculations. Within the approximations made the agreement is 
encouraging. 



4 ROTOVIB RATIONAL SPECTRUM 

Given the aligned and transverse vibrational energies, we construct 
the rotovibrational spectra of helium molecules. From the equations 
lf5J and ([6), the molecular system has a finite zero-point energy 
associated with aligned and transverse vibration 



f000 = ^(hcj, + hCQu) 



(25) 



when N = A = V = 0. Therefore, the actual dissociation energy 
will be reduced from D,„ by e oo (LS96). Figure [3] shows the rotovi- 
brational energy spectrum of He2 in the ground state at B = 10 12 
G and 10 13 G respectively. Table [8] shows the number of rotovi- 
brational states of various helium molecular ions. The number of 
rotovibrational states decreases for higher electronic excited states. 
The magnetic field dependence is more complicated for the fol- 
lowing reasons. The number of aligned vibrational levels generally 
increases with B. On the other hand, the number of transverse vi- 
brational (A 7 ) and rotational (A) levels increases with B at B S Bq 
(B q = 4.414X 10' 3 G) while it decreases with B at B>B Q . This is be- 
cause the ion cyclotron energy (oc B) dominates over ha> ± (oc B >/2 ) 
at B^ Bq. It should be noted that only those rotovibrational states 
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Table 8. Number of rotovibrational states in helium molecular ions. 
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Figure 3. Rotovibrational energy spectrum (ejvAv) of He2 in the ground 
state at B = 10 12 G (top) and 10 13 G (bottom). The zero energy corresponds 
to E = —E m , while the uppermost horizontal line corresponds to E = —E m + 
D m above which the molecular states become dissociated. The black, blue 
and red lines indicate N = 0, 1 and 2 state. For each N and A, the horizontal 
lines indicate energy levels for different V states. 



with excitation energies that are smaller than or similar to the ther- 
mal energy have large statistical weights in the partition function 



5 EFFECTS OF FINITE NUCLEAR MASS 

The separation of the cent er-of-mass motion is non-trivial when 
a magnetic field is present jHeroldetal]|l98ll) . The total pseudo- 
momentum (K) is often used to take into account m otional eff ects 
in a magnetic field since a is a constant of motion 01200 lh . In 
the following subsections, we discuss two effects associated with 
finite nuclear mass in strong magnetic fields. We denote the bind- 
ing energy in the assumption of fixed nuclear location (e.g., Born- 
Oppenheimer approximation) as ef\< 0) for an electronic state 
k. Since we consider states with n = and v = 0, the relevant 
quantum numbers are k = (m,j (( denotes each bound electron in 
multi-electron atoms and molecules). 



5.1 Finite nuclear mass correction 

The assumption of zero transverse pseudomom entum introduces an 
additional term s,;hQ. B in the binding energy dHerold et al.ll 198 ll ; 



IWunner et al.ll 198 Th . Q.b is the nuclear cyclotron energy and s K = 
2i >Wj is the sum of magnetic quantum numbers for a given elec- 
tronic state k (e.g., s K = 4 for He 2 molecule in the ground state). 
However, the scheme assuming the zero transverse pseudomo- 
mentum does not necessarily give the lowest binding energies at 
B> Bq. Instead, LS95 and LS96 estimated lower binding energies 
at B 2; Bq using another scheme which relaxed the assumption of 
the zero transverse pseudomoment um. A more rigorous cal cula- 
tion was performed for He + i on by iBezchastnov et al.l i ll 9981) and 
IPavlov & Bezchastnovl J2005b . An application of such schemes to 
multi-electron systems is beyond the scope of this paper. We will 
discuss the limitation of our models at very high magnetic field in 



5.2 Motional Stark effects 

When an atom or molecule moves across the magnetic field B, a 
motional Stark electric field E MS = is induced in the center-of- 
mass frame. M is the mass of atom or molecule. The Hamiltonian 



for the motional Stark field is given by 



Hi, 



KxB 
' ; 

Mc 



(26) 



where fi s0 = Jjp. For a given pseudomomentum ft, the motional 
Stark field separates the guiding center of the nucleus and that of 
the electron by 



Rk 



c\K x B\ 
eB 2 



cK ± 
eB ' 



(27) 



Since the motional Stark field breaks the cylindrical sym- 
metry preserved in magnetic field, it is non-trivial to evaluate 
motional Stark field effects. A non-perturbative (therefore more 
rigoro us) approach has been applied only for one-electron sys- 
tems dVincke et all 1 19921 ; IPotekhinl 1 1994 IPavlov & Bezchastnovl 
120051) . However, such an approach is quite complicated and time- 
consuming especially for multi-electron atoms and molecules. 
Therefore, following LS95, we considered two limiting cases: (1) 
R K <k p and (2) R K » p and determined general formula which can 
be applied to a wide range of B and K ± . For diatomic molecules, we 
can apply a nearly identical scheme used for calculating transverse 
vibrational energy to the both cases. 
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5.2.1 Centered states 

When the energy shift caused by motional Stark field is smaller 
than the spacing bet ween binding energ i es, th e perturbation ap- 
proach is applicable l lPavlov & Meszaroa j 1993b and § 13. 2. It , The 
first order perturbation energy (k\H M s \k) vanishes since the matrix 
element (k\x\k) = 0. The second order perturbation energy is given 
by 



J2) 



z 



\W\H ms \k)\ 2 
JO) _ (0) 



K l (2) 
u 

2M * 



(28) 



Among the various electronic states k' = {n'm'v'}, only n' = 0, m' = 



m± 1, V = state have non-negligible contribution to e^ 2 ' since the 



other states have large e^? - e"" and/or vanishing matrix element 
{k'\x\k). Since the overlap integral of lon gitudinal wavefunction is 
close to unity dPavlov & Mes zaros 1993]), o^ 2 ' is given by 



JO) 



v(2> 



hn B 



m + 1 



' £m+l + hQ-Bi fim _ £»t-] _ hQ.B 



(29) 



where e„, = -ef\> 0) where m denotes the magnetic quantum 
number of the outermost electron of ele ctronic state k. The 2nd term 
is zero for the ground state. Following IPavlov & Meszarosl Jl993h . 
we define the anisotropic mass M ±JI as 



Mi. = 



M 



1 



> M. 



The transverse energy characterized by K ± is given by 

K ± 2 



2M ± / 



(30) 



(31) 



The perturbation m ethod is valid when \e[ 2) \ <k |Ae^ 0) | 



dPavlov & Meszarosf l 993 ) where K ± K is given by 

/2M|A4°M 1/2 



(32) 



where |As£ | is the spacing of the zeroth order energies (typically 



5.2.2 Decentered states 

When R K » p, it is convenient to utilize the so-called decentered 
formalism (LS95). We replace the nuclear Coulomb term by 



R K ) = J d 2 ? ± 



l*„,(^)l 2 



(33) 



\?±+Rk\ 

and compute binding energies at different R K grid points. For di- 
atomic molecules, we replace the nuclear Coulomb term by 



(34) 



V(z,R K ) = V m (z-^,R K ) + V m (z + ^,R K ). 

The grid calculation for R K is identical to the one for transverse vi- 
brational energy in j j3.2| except that the Coulomb repulsion term be- 
tween two nuclei is (instead of t Z T ,„ ) since motional Stark 

"0 v (05 +R 2 K )" 2 ' 

field shifts the guiding center of the two nuclei by R K in the trans- 
verse direction but the separation between the two nuclei is still 

a . 

LS95 found that the binding energy curves are well fit by the 
following formula. 



E ± (K ± ) + e 



(0) 



-At In 



1 



A 2 + A 3 Rl 



(35) 



> 

CD 








Figure 4. The transverse energy of He2 molecule at B = 10 G. The as- 
terisks are the binding energies from our numerical calculation and they are 
fitted with the function given in < I35 1 . The dashed line shows the energy 
curve from the perturbation method discussed in j|3.2. H and 35.2.11 



where A\,Ai and A3 are the fit parameters. E ± (R K ) = .EjXA'j.) is the 
transverse energy and e^ 0> is the binding energy in the infinite nu- 
clear mass approximation. Figure|4]shows the binding energy curve 
of He 2 molecule as a function of R K at B = 10 12 G. At small R K , 
the fitted function is well matched with the results from the per- 
turbation approach. Althoug h mixing betwee n dif ferent m s t ates is 
ignored, a comparison with lPotekhinl dl998h and IPotekhir] ( 1 1994 
for hydrogen atoms indicates this a pproach gives better t han 30% 
accuracy over a large range of K L dLai & Salpetej[T995l) . This is 
adequate for our purpose of investigating ionization and dissocia- 
tion balance. 

Along with the finite nuclear mass term discussed in 35.11 the 
electronic energy of an atom or molecule moving with transverse 
pseudomomentum K ± is given by 



e,(#x) = 4°' + s K h£l B + Ej_(K x ). 



(36) 



Note that e K < and both the second and third term decrease the 
binding energy. Later, we will discuss the validity of this approach 
at B>B Q . 



6 IONIZATION AND DISSOCIATION EQUILIBRIUM 

We have investigated the ionization and dissociation balance of 
magnetized helium atmospheres including the following chemical 
reaction channels. Table |9l and 1101 list ionization and dissociation 
energies of various helium ions and molecular ions in the assump- 
tion with fixed nuclear location. We did not take into account the 
He" ion since its ionization energy is < 20 [eV] and He" is not 
abundant at all in the temperature range considered here (T > 10 5 
K). 

• Ionization 

(1) He + <-> a + e 

(2) He ^ He + + e 

(3) He^ ^ He' + + e 

(4) He+ «^He 2+ + e 

(5) Hej ^He+ + e 

• Dissociation 

(6) He^ + «->q' + He + 
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Table 9. Ionization energy [eV] of helium atom and molecular ions. 



Bl 2 


He+ 


He 


He 3+ 


He 2+ 


He 2 


He 2 


1 


418 


159 


391 


418 


261 


137 


10 


847 


331 


885 


947 


603 


298 


100 


1563 


629 


1833 


1912 


1209 


643 



Table 10. Dissociation energy D„, [eV] of helium molecular ions. The num- 
bers in the brackets are D,„ [eV]. The columns without a number indicate 
that there is no bound molecular state with respect to ions and atoms. The 
zero-point energy correction is not included. 



Bi2 


He 3+ 


He 2+ 


He 2 


He 2 


1 


-(-) 


-(97.8) 


74.6(251) 


53.1 (290) 


10 


37.8 (37.8) 


137 (359) 


410 (725) 


378 (806) 


LOO 


270 (270) 


619 (973) 


1198(1708) 


1212(1915) 



Notes: Recent work by Turbiner & Lopez Vievral J2004T) and 
iTurbiner & Guevari] )2006h shows He 3+ and He 2+ ion can exist at 
B > (2.4 - 2.6) X 10 12 G. Our results are roughly consistent with theirs. 



(7) He 2+ <-» a + He 

(8) He 2+ <-> He + + He + 

(9) HeJ <-> He + + He 
(10) He 2 <-> He + He 



The Saha-Boltzmann equation for ionization balance is given 



as. 



n M n e Zi+iZe 

where z e is the partition function for an electron, 



' m e kT 
i 2nh 2 



3/2 



tanh T] e 



(37) 



(38) 



r] e = h u> B /2kT and a> B = eB/m e c is the electron cyclotron fre- 
quency dGnedin et alj|l974h . The quantity z, is the partition func- 
tion for ionization state i, 



' MjkT 
, 2nh 2 



1/2 



where Ej, Z ; and Mj are the ground state energy, the charge and the 
mass of an ion i. 77, = hSl Bj /2kT where Q. Bj = ZjeB/MjC is the ion 
cyclotron frequency. ei, K (K ± )(< 0) is the binding energy of a bound 
state k given in equation J36b . w,> is the occupation probability. In 
general, w, A ,(A* ± ) is a function of Z,, Z p (the charge of perturbing 
ions, usually the effective charge of the plasma), ej K (K L ) and p (the 
plasma density), w depends not only on the type of ion (', the elec- 
tronic state k and the transverse pseudomomentum K±, but also it is 
different for each bound electron. Obviously, electrons in the outer 
shells are subject to a stronger electric field from neighboring ions 
than those in the inner shells. We explicitly computed w, A for all 
the bound electrons u sing the electronic microfield distribution of 
IPotekhin et all d2002l) . 

The Saha-Boltzmann equation for dissociation balance is 
given as 



n k n, 



ZkZl 



(40) 



where Zj is the partition function for molecular ionization state j, 
'MjkT) 



12 



2nh 2 } sinh 



1j 



1 mfZ^ 



x exp 



AT 



e NAV <D k 



2j ex P(- 



kT I 



(41) 



where Ej, Z, and Mj are the ground state energy, the charge and the 
mass of a molecular ion j. rjj = hQ. B j/2kT where Q B j = ZjeB/MjC 
is the molecular ion cyclotron frequency. e,>(K ± )(< 0) and Wj K (K ± ) 
are the binding energy and occupation probability of an electronic 
bound state k. £nam(> 0) is the excitation energy of a rotovibra- 
tional state (N, A, V). We took the summation zZna v unt i' e NAV ex- 
ceeds the dissociation energy D K . The set of the Saha-Boltzmann 
equations along with the condition for the baryon number conser- 
vation and charge neutrality are iteratively solved until we reached 
sufficient convergence in An e /n e (< 10 -6 ) where n e is the density of 
free electrons. The convergence was achieved rapidly in most cases 
(less than 10 iterations were required). 

Figure|5]|6]and|7]show the fraction of helium ions and molecu- 
lar ions at B = 10 12 , 10 13 and 10 14 G. At B = 10 12 G, He 3+ and He 2+ 
are not present because they are not bound with respect to their dis- 
sociated atoms and ions (table [Tot. In all the cases, the He?, + frac- 
tion is negligible because helium molecular ions with more elec- 
trons have much larger binding energies. At B = 10 14 G, He 2 is 
not present because even the ground state becomes auto-ionized to 
HeJ ion due to the finite nuclear mass effect although He 2 remains 
bound with respect to two helium atoms. Note that the ionization 
energy of He 2 at 10 14 G is 643 eV (table|9j while the difference in 
the energy due to the finite nuclear mass term between He 2 and HeJ 
is 945 eV. As discussed in 35.11 our scheme using the zero pseudo- 
momentum becomes invalid at B £ Bq and in reality He 2 will have 
lower binding energy than HeJ . Therefore, our results at B = 10' 4 
G should be taken with caution. Figure [8] shows the temperature 
dependence of the helium atomic and molecular fractions at dif- 
ferent B-field strengths. We fixed the plasma density to the typical 
density of the X-mode phot osphere (~ 10, 10 2 an d 10 3 g/cm 3 for 
B = 10 12 , 10 13 and 10 14 G jLai & Salpeterll 19971) ). The transition 
from molecules to helium atoms and ions takes place rapidly at 
T ~ 3 x 10 5 K and ~ 6 x 10 s K for B = 10 12 G and 10 13 G. 

In order to illustrate different physical effects, we show the 
fraction of helium atoms and He 2 molecules as a function of tem- 
perature in Figure|9] When motional Stark effects are ignored (dot- 
ted line), the molecular fraction is underestimated because helium 
ions and atoms are subject to a larger motional Stark field due to 
their smaller masses and binding energies. As expected, molecules 
become more abundant when rotovibrational states are included. 
The discrepancy from the results without rotovibrational states 
(dashed line) increases toward higher temperature since more ro- 
tovibrational states have larger statistical weight as their excitation 
energies are an order of 100 eV at B = 10 13 G (see figure|3j. When 
we take into account only the ground states (dot-dashed line), the 
results are close to the case without motional Stark effects because 
the ground states are less affected by motional Stark effects than the 
excited states. 



6.1 Hes molecule and larger helium molecular chains 

When He 2 is abundant, larger helium molecules such as He 3 
may become abundant. We roughly estimated the fraction of He 3 
molecules by neglecting the finite nuclear mass effects and zero- 
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Figure 5. Ionization and dissociation balance of helium at T = 10 5 5 K (top) 
and 10 6 K (bottom) at B = 10 12 G. 
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Figure 6. Ionization and dissociation balance of helium at T = 10 5 - 5 K ( top) 
and T = 10 6 K (bottom) at B = 10 13 G. 



point vibrational energy correction both of which will reduce the 
dissociation energy. Similarly to He2, we computed the binding en- 
ergy of He3 molecules by Hartree-Fock calculation. The dissoci- 
ation energy of He 3 ^He 2 +He is 289 and 1458 eV at B = 10 13 
and 10 14 G. Our results are close to those of the density functional 
calculation by ML06 (384 and 1647 [eV] at B = 10 13 and 10 14 
G). Note that the density function calculation overestimates bind- 
ing energies by ~ 10% compared with more accurate Hartree-Fock 
calculation (ML06). Assuming the internal degrees of freedom are 
similar between He 2 and He3 molecule and neglecting the bending 
degree of freedom for He 3 , we computed dissociation equilibrium 
between He2 and He3 following lLai & Salpeted i 19971) . 

Figure [Tol shows contours of the He3 molecule fraction with 
respect to He 2 molecule at B = 10 13 G. In the dotted area, the 
fraction of diatomic helium molecular ions is larger than 10%. It is 
seen that He3 molecule will be abundant at T < 3 x 10 5 K at B = 
10 13 G. The result s are consis tent with the estimates of molecular 
chain formation by [Eaild200lh and ML06. ML06 provided binding 
energy of infinite He molecular chain as E m a 1 .256j 3 38 keV at B > 
10 13 G. From this fitting formula, the cohesive energy of helium 
molecular chains is given as E co ~ 0.36, 1.5 and 5.1 keV at B = 
10 13 , 10 14 and 10 15 G. W hen k T < 0AE co , molecular chains are 
likely to be formed JEaill200lh : consequently, we expect helium 
molecular chains to form at T < 3 x 10 5 , 1 x 10 6 and 5 x 10 6 K and 
B > 10 13 , 10 14 and 10 15 G. 



7 APPLICATION 

We investigated the ionization and dissociation balance of helium 
atmospheres for two classes of INS whose X-ray spectra show ab- 
sorption features. 

7.1 Radio-quiet neutron stars 

A class of INS called radio-quiet neutron stars (RQNS) is char- 
acterized by their X-ray thermal spectra with T < 10 s K and 
spin-down dipole B-field strength B > 10 13 G. A single or mul- 
tiple jtbsor2tiojifeatoesJ^e bera_^ radio-quiet 
NS dHaberl et alj|2003ll2004|jvan Kerkwiik et al.ll2004l) . Although 
the interpretation of these features is still in debate, helium is cer- 
tainly one of the candidates for the surface composition of RQNS 
dvan Kerkwiik & Kaplanl |2006|) . Timing analysis suggests that a 
couple of RQNS ha ve dipole magnetic field streng ths in the range 
of B = 10 13 -10 14 G dKaplan & van Kerkwiikll2005l) . Therefore, we 
investigated the ionization/dissociation balance of a helium atmo- 
sphere at B = 3 x 10 13 G. Figure QT] shows the contours of He + and 
He molecular fractions at 3 x 10 13 G. He + is dominant at T ~ 10 6 
K and molecules become largely populated at T < 5 x 10 5 K. 

Suppose all the RQNS have helium atmospheres on the 
surface. RQNS with higher temperatures (T ~ 10 6 K) such 
as RXJ0720.4-3125 and RXJ1605.3+3249 will have He + ions 
predominantly with a small fraction of He atoms. Indeed, sev- 
eral bound-bound transition lines of He + ion have energies 
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Figure 7. Ionization and dissociation balance of helium at T = 10 5 5 K (top) 
and T = 10 6 K (bottom) at B = 10 14 G. 



at the observed absorption line location at B > 10 13 G 
dPavlov & Bezchastnovl 120051) . On the other hand, RQNS with 
lower temperatures (T~5x 10 5 K) such as RXJ0420.0-5022 may 
have some fraction of helium molecules, although the He + ion is 
still predominant at low densities where absorption lines are likely 
formed. 



= 10 log(p) = 2 B 12 =100 Iog0o) = 3 
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Figure 8. Temperature dependence of helium atomic and molecular fraction 
at B = 10 12 G (p = 10 g/cm 3 ), B = 10 13 G (p = 10 2 g/cm 3 ) and B = 10 14 
G (p = 10 3 g/cm 3 ). 
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Figure 9. Fraction of helium atoms (left) and helium molecules (right) 
p = 10 2 g/cm 3 at B = 10 13 G. The solid lines show the results by tak- 
ing into account all the physical effects discussed in this paper. The dotted 
and dashed lines ignored motional Stark effects and rotovibrational states 
respectively. The dot-dashed lines include only the ground states. 



7.2 1E1207.4-5209 

1E1207.4-5209 is a hot isolated NS0 with age ~ 7 x 10 3 yrs. Th e 
fitted blackbody temperature is ~ 2 x 10 6 K jMori et alj|2005l) . 
Presence of a non-hydrogenic atmosphere has been suggested 
since 1E1207 shows multiple absorption featur es at higher energies 
than hydrogen atmo s phere models predicte d (Sanwal et al J 1 20021; 
Hailev & Morill2002| ; IPavlov & Bezchastnovll200^; iMori & Haiievl 
20061) . ISanwal et all J2002h interpreted the observed features as 
boun d-bound t r ansitio n lines of a He + ion at 6 = 2 x 10 14 G. 
Also, iTurbinej d2005f) suggested that He 3+ molecular ion may be 
responsible for one of the absorption features observed in 1E1207 
atB~4.4x 10 13 G. 

However, our study shows that the fraction of He 3 / molecu- 
lar ions is negligible at any B-field and temperature because more 



2 Recent timing a nalysis suggests that 1E1207 is in a binary system 
iWoods et ai]|20O6h . However, it is unlikely that 1E1207 is an accreting 
NS. 



neutral molecular ions have significantly larger binding energies. 
We also investigated ionization/dissociation balance of helium at- 
mospheres at 2 x 10 14 G (figure [T2l. At B <T Bq, our scheme of 
treating finite nuclear mass effects becomes progressively inaccu- 
rate. As a result, the ionization energy of helium atoms becomes 
significantly smaller and the He2 molecule becomes auto-ionized 
due to the finite nuclear mass correction discussed in 35.11 There- 
fore, the molecular fraction will be underestimated in this case (the 
left panel in figure \\7\. On the other hand, when we ignore the 
finite nuclear mass effects (the right panel in figure [T2l. the molec- 
ular fraction is likely overestimated. We expect that a realistic ionic 
and molecular fraction will be somewhere between the two cases. 
It is premature to conclude whether He + is dominant for the case 
of IE 1207 since we do not have a self-c onsistent te mperature and 
density profile. The study of ML06 and lLail d200lh also suggests 
the critical temperature below which helium molecule chains form 
is ~ 2 x 10 6 K at B = 2 x 10 14 G. This is close to the blackbody 
temperature of 1E1207. Further detailed studies are necessary to 
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Figure 10. Ratio of He3 with respect to He2 at B = 10 13 G (solid curve con- 
tours). The dotted area shows the (p, T) regime where a fraction of diatomic 
helium molecular ions is larger than 10%. 




log T [K] 

Figure 11. Ionization and dissociation balance of helium at B = 3 X 10 13 
G. The dotted area is where the fraction of helium molecules and molec- 
ular ions (including He2, HeJ, He 2+ and He 3+ ) becomes more than 50%. 
The black area is where the fraction of He + becomes more than 50%. The 
white region on the right is where the fraction of bare helium ions becomes 
more than 50%. In the narrow white region betweeen the dotted and black 
area, atomic helium is somewhat abundant, so neither molecular nor singly 
ionized helium dominates. 



conclude the composition of helium atmospheres at B = 2 x 10 
G. 



8 DISCUSSION 

We have examined the ionization-dissociation balance of the 
helium atmospheres of strongly magnetized neutron stars. As 
the observational data on isolated neutron stars has improved 
over the past decade {Sanwal et alj 120021 : lHailev & lyiorfl 20021 ; 
Ivan Kerkwiik & Kaplarj l2006h . both hydrogen and iron atmo- 
spheres have lost favour, and atmospheres composed of other ele- 
ments have been invoked to interpret both the continuum and spec- 
tral properties o f neutron sta r s. On the other hand the theoretica l 
investigations of IChang etal] d2004h and lChang & Bildstenl | |2004|) 
have shown that diffusive nuclear burning may quickly deplete the 




log T [K] 

Figure 12. Ionization and dissociation balance of helium at B = 2 X 10 14 
G. The top and bottom figure show the results with and without the finite 
nuclear mass effects respectively. In the dotted region the molecular and 
molecular ion fraction exceeds one half. In the black region singly ionized 
atomic helium dominates, and the white region demarcates fully ionized 
helium. 

hydrogen by so the NS surface may be composed of helium or 
heavier elements. 

The results found here indicate several avenues for future 
work. Although the treatment of finite nuclear mass effects is prob- 
lematic in the strong field limit, only by including this accurately 
can we know the physical state of matter in super-critical magnetic 
fields - solid, liquid or gas. The state can have significant effects 
of the outgoing radiation from these objects. It is also important 
to repeat a similar calculation to that presented here for molec- 
ular chains, including all of the relevant degrees of freedom of 
the chains, i.e. including the bending modes which were neglected 
here. Regardless, the high abundance of molecules under the con- 
ditions of observed neutron star atmospheres is bound to spur addi- 
tional research into the statistical mechanics of highly magnetized 
molecules and polymers. 

At the temperatures and densities of neutron star atmospheres 
the rotovibrational excitations of helium molecules are populated. 
Including these excitations increases the expected abundance of 
molecules by up to two orders of magnitude relative to calculations 
that ignore the internal states of the molecule. Ionization, dissocia- 
tion and electric excitation energies of helium molecules are larger 
than 100 eV at BH 10 13 G. On the other hand, rotovibrational exci- 
tation energies are in the range of 10-100 eV at B = 10 12 — 10 14 G. 
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If helium molecules are abundannt, their spectroscopic signatures 
may be detected in the optical, UV and X-ray band. If helium com- 
prises the atmospheres of isolated neutron stars, clearly it is cru- 
cial to understand the structure of helium molecules and molecular 
chains in order to interpret the spectra from neutron stars. 
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